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O ■ ABSTRACT 
O ■ 

<N : 

, We propose a method for estimating the large-scale rms bulk velocity of the 

O ! cosmic mass field from the transmitted fluxes of Lya forests. It is based on two 

lO ! linear relationships on large scales: 1) the relation between the fluctuations of the 

^ I transmission and the underlying density field, and 2) the relation between the 

^ I density fluctuations and the peculiar velocity field. We show that, with a multi- 

i/-^ [ scale decomposition, the two relations can be effectively employed for predicting 

^ I the rms bulk velocity. Since QSO's Lya forest is due to the absorptions of dif- 

^ I fusely distributed and photoionized IGM, this method provides an independent 

^ I estimate of the rms bulk velocity at high redshifts, on large scales, and free from 

^ I the bias of galaxies. Using the transmitted flux of 60 moderate-resolution QSO 

Ph| spectra, the rms bulk velocity is found to be 230±50 km s~^ around redshift 

O I z = 2.25 on scale 23 h~^Mpc, and down to 110±45 km s~^ around z = 3.25 on 

I scale 92 h~^Mpc for an LCDM universe {Q = 0.3 and A = 0.7). The results are 

' basically consistent with the linear evolution theory. 

> : 

^ I Subject headings: cosmology: theory - large-scale structure of universe 



1. Introduction 



The transmitted fluxes of high redshift QSOs' Lya absorption spectra have been success- 
fully applied to many aspects of cosmic large-scale structure study, such as the discrimination 
among dark matter models (Bi, Ge & Fang 1995; Bi & Davidsen 1997), recovery of the initial 
linear mass power spectrum and estimation of cosmological parameters (Croft et al. 1998 & 
1999; Hui 1999; Nusser & Haehnelt 1999; McDonald & Miralda-Escude 1999; Feng & Fang 
2000; McDonald et al 2000), and finding the applicable range of the hierarchical clustering 
model (Feng, Pando & Fang 2001; Zhan, Jamkhedkar & Fang 2001). These studies show 
that Lya forests provide an important complement to structure formation studies based on 
galaxy samples. 
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In this paper we extend the apphcation of Lya forests to the cosmic velocity field. 
In terms of cosmological parameter determination, the rms bulk velocity is an important 
statistic, as it is related to an integrated power spectrum over the scales beyond the directly 
detectable range of the power spectrum with current galaxy surveys. On such scales, the 
cosmic clustering remains in the linear regime, and, therefore, would be useful to set a 
constraint on cosmological parameters, and on the relation between the galaxy distribution 
and the underlying mass field. 

Many efforts have been made on measuring and estimating cosmic bulk velocity with 
galaxy samples. An important result is that the velocity field calculated from spatial distri- 
bution of galaxies is found to be basically consistent with the peculiar velocities measured 
based on the TuUy-Fisher or fundamental plane relations (e.g. Dekel et al 1999; Branchini 
et al 1999). This indicates that the bulk velocity of galaxies on large scale is of gravitational 
origin. That is, the cosmic velocity field v(x) on large scales is related to the cosmic mass 
density field 5(x) via the following hnear relation or its variants 



where / ~ for the local universe. We define x = {x,y,r), where x and y are the 
coordinates on the celestial sphere, and r is the radial or redshift direction. 

A problem with galaxy samples, however, arises from the likely bias of the distribu- 
tions of galaxies relative to the underlying mass distribution in the universe. The simplest 
bias model assumes that the number density fluctuations of galaxies, Sg{x), is related to 
the underlying mass density field by 5g(x) = 65 (x), where b is the bias parameter. More 
sophisticated models assume that the relation between ^^(x) and (5(x) might be stochastic, 
non-local, and non-linear. Furthermore, the peculiar velocities of galaxies may also be biased 
from the cosmic velocity field. These effects will cause uncertainty in deriving the cosmic ve- 
locity field from galaxy samples. It is therefore important to develop a method of estimating 
the cosmic velocity field free from galaxy bias. 

Here we calculate, based on Eq. (1), the large-scale rms bulk velocities from samples 
of the transmitted flux F of QSO's Lya absorption spectrum. The QSO Lya forest is 
very well modeled by absorption from diffusely distributed and photoionized IGM, which 
in turn is believed to be distributed proportionally to the density of the underlying mass 
field p(x) on scales larger than the Jeans length of the IGM (for a review, see Rauch 1998). 
The transmitted flux of QSO's Lya forest is thus likely to be less biased as a tracer of 
the underlying mass and velocity field, providing an independent estimate of the rms bulk 
velocity, separate from that defined by galaxies. 

To do so, we need to develop methods for: 1) Obtaining mass fluctuation data 5(x) from 
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Lya flux F in QSO spectra, in the linear or quasi-linear regimes. 2) Properly calculating 
the bulk velocity from mass fluctuations. These two methods are discussed, respectively, in 
§2 and §3. The results of our investigation of the rms bulk velocity on large scales at high 
redshifts are presented in §4. 



2. Probability Distribution Function of Lya Transmission 

2.1. Transmission— density relation 

For a QSO with coordinate on the celestial sphere, the observed flux of Lja 

absorptions is equal to Fce~'^^^\ with Fc being the continuum, and t{z) the optical depth. 
The data to be used have been normalized, and therefore the transmission F(z) = e~'^^^\ 

For a diffusely distributed HI cloud in photoionization equilibrium, we have (e.g. Bi 
1993; Fang et al. 1993) 

t{z) = A{z) J dr[l + 5{xi, Vi, r)]''V[w{z) - r - Vrixi, yi, r), b] (2) 

where w{z) is the redshift space coordinate, r is the radial comoving coordinate, V is the 
normalized Voigt profile that is approximately Gaussian, V ~ l/{y/7rb) x exp{— [^(2;) — r — 
■^^^(2^1; Z/i) 0]^/^^}) '^rixi-iUi-iT) is the radial component of peculiar velocity, h is the Doppler 
broadening on the order of several tens km s~^, and 5{xi, i/i, r) = [p{xi, yi, z) —p]/p is the dark 
matter density contrast smoothed on the scale of the IGM Jeans length (e.g. Bi & Davidsen 
1997; Nusser & Haehnelt 1999). The parameter A{z), which depends on the cosmic baryonic 
density, the photoionization rate of HI, and the mean temperature of IGM, is referred to 
the mean transmission over a redshift range as {F{z)) = (e""^*^^^) — e~^^^). The parameter a 
depends on reionization, and is in the range of 1.5 to 1.9 (Hui & Gnedin 1997). 

There arc two problems in drawing information of linear or quasi-linear density fluctua- 
tions 6 from the transmitted flux F. First, the relation between F and 6 is non-linear. This 
is especially true for a saturated absorption hne, i.e. F ~ corresponds to the non-hnear 
regime of 5. Second, the peculiar velocity Vr, which is what we are trying to detect, already 
enters in Eq.(2). These two problems can be avoided by considering the smoothed F and 
5 on large scales. If the smoothing scales are much larger than the Doppler factor b and 
the positional uncertainty caused by the peculiar velocity v{x), the Voigt function V can be 
approximated as a Dirac delta function d^{r — w{z)). In this case, Eq.(2) yields 



F{z)c:^exp{-A{z)[l + 5{x„y„z)r}. 



(3) 
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For simplicity, in Eq.(3) and hereafter we still use the notations F{z) and 5{z) for the 
smoothed transmission and density perturbations. Thus, for the linear or quasi-linear regime, 
i.e., \SF{z)\ = \F{z) - {F{z))\/{F{z)) < 1 and \5{xi,yi,z)\ < 1, Eq.(3) gives 

6{xi, Vi, z) ~ —-^5F{z), (4) 

where {xi,yi) is to remind that it is for the i-th QSO. 

Eq.(4) shows that the fluctuations of the Lya transmission flux 5F{z) trace the underly- 
ing density perturbation 6{xi, yi, z) point-by-point. The factor aA{z) plays the role of a bias 
parameter between the Lya flux fluctuations and the mass density perturbations. aA[z) is of 
order of 1, as a ~ 1.5 — 1.9 and A{z) ~ 0.75 — 0.9. Eq.(3) has been used in the recovery of the 
initial linear power spectrum of the underlying mass field by the observed Lya transmission 
(e.g. Croft et al. 1998). It should be reminded that Eq.(4) holds only if 1) the variance of 
5F is less than 1, or the variance of F is less than {F{z)); and 2) the positional uncertainty 
caused by the peculiar velocity is much smaller than the scales considered. We will check 
the validity of these conditions in the analysis below. 



2.2. Transmission— density relation in DWT representation 

To do the smoothing of the transmitted fluxes, we apply the algorithm of multiscale- 
decomposition based on the discrete wavelet transform (DWT) (Daubechies 1992; Fang & 
Thews 1998). Because the bases of the DWT analysis are orthogonal, complete and localized, 
it is easy for the calculation to use Eq.(l) (next section). In our calculation, we use the 
Daubechies 4 (D4) wavelets. 

Let us consider a sample 5F{z) that spans a redshift space L = Z2~zi, and the spectrum 
is binned into 2'^ pixels with J being an integer. With the DWT analysis, the smoothed 
SF{z) on the scale is given by 

2^-1 

5Fjiz) = J2ef,<Pj,{z), (5) 

1=0 

where (f)ji{z)s are the DWT scaling functions for mode (jl). The label j is for spatial scale 
, and / is for the position in redshift space around zi + IL/2K The scaling function plays 
the role of a window function. The scaling function coefficient (SEC) ej^ is given by 



4ii^i^yi)= / SF{z)(pji{z)dz, 



(6) 
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where {xi, yi) is to show that the flux fluctuation 5F{z) in Eq.(6) is for the i-th QSO. The 
SFC e^i is proportional to the mean flux over a bin size around the position I. Thus, 
with the DWT decomposition, Eq.(4) becomes 

where the SFC of the mass fleld is given by 

^Tii^i^ Vi) = j z)(j)ji{z)dz. (8) 

Therefore, Eq.(7) is the DWT representation of Eq.(4). The SFCs e|J and actually are, 
respectively, the fluctuations SF{z) and S{z) in DWT representation. Eq.(7) has been used 
in the DWT recovery of the initial linear power spectrum of the underlying mass fleld by the 
observed Lya transmission flux (Feng & Fang 2000). 



2.3. PDF of Lya Transmission 

The samples of transmitted fluxes F{z) used in our study comprise of 60 QSOs' spectra 
selected from Bcchtold (1994), Dobrzycki & Bechtold (1996), Scott, Bechtold & Dobrzycki 
(2000), and Scott ct al. (2000). The selection ensures that each spectrum in the calculation 
below is comparable to at least twice the length of the largest scale (92 h~^Mpc) interested. 
The emission redshift of the QSOs covers the range from 1.9 to 4.12. Each spectrum is 
averaged in bins of size ~ 1.2 A, such that there are 1024 'pixels' per unit redshift, which 
converts to a resolution of c/1024 ~ 293 km s"^. 

Using Eq.(6), we compute the SFCs ej^ of flux fluctuations 6F{z) for each QSO. We use 
the so-called off-counting method to treat the spectra (Jamkhedkar, Bi & Fang 2001), i.e. 
the SFC e^i of modes (j, /) is not counted in the calculation if the data at that mode are 
contaminated by metal line, bad, or low S/N pixels. An advantage of the DWT off-counting 
method is that the statistical results on large scales are not sensitive to the bad pieces on 
small scales. 

We treat all the data as ensembles of various modes (j, /). That is, for a given mode 
{j, I), the ensemble consists of all SFCs ej^(xi, yi) of the QSOs that have flux observed around 
A = (1 -I- 2^1 -I- lL/2^)Xa, where ~ 1216A. Thus, we can calculate the ensemble average of 
SFCs of each mode (j, I) by 

1 ^ 

i=l 
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where N is the number of the SFCs in the ensemble. 

Since modes (j, /) are locahzed in redshift space, we can also construct ensembles of 
the SFCs in different redshift ranges. For a given redshift range z ± {Az/2), the ensemble 
consists of all eji{xi, yi) of which the label I corresponds to position in the range z ± {l^z/2). 
In the following calculations, we use three redshift ensembles. They are z—2.2b, 2.75, and 
3.25 with = 0.50. 

The probabihty distribution functions (PDFs) of the SFCs ej^ are calculated in the three 
redshift ranges. The result for z — 2.75 is plotted in Figure 1. The scale j corresponds to 
redshift distance 5z = Az/2^ — 2~^~^. In Fig. 1, the physical scale of j is calculated from 
the redshift distance 5z z — 2.75, for a LCDM model (Q = 0.3 and A = 0.7). For instance, 
the scale 92 h~^ Mpc corresponds to 5z = 0.125, or in wavelength AA ~ 152A. 

In Figure 1, the panel on the smallest scale, i.e. 72 h~^kpc, is the PDF of the binned 
distribution F{z) without the SFC smoothing, and therefore, it actually is the PDF of QSO 
Lya transmission, which has been done by many others (e.g. McDonald et al. 2000). Rather, 
Figure 1 shows the scale dependence of the PDFs. When the smoothing scale is small (large 
j), the PDF is significantly non-Gaussian, while it is essentially Gaussian on large scales. 
The K-S test shows that the PDFs are Gaussian on scales larger than 12 h~^ Mpc, and non- 
Gaussian on scales less than ~ 12 h~^ Mpc. This is consistent with the authors' previous work 
(Zhan, Jamkhedkar, & Fang 2001). Moreover, the variance of 2(-^~^)/^ej^, or equivalently the 
variance of 5Fj{z), is less than 1. Thus, it is reasonable to calculate the density perturbations 
of the underlying mass field by Eqs.(4) or (7). 

All the above-mentioned properties are typical. That is, in other redshift ranges, the 
behavior of the ej^ PDFs is similar to Figure 1 with a slight evolution in the redshift range 
considered. Thus, on scales larger than 12 h~^ Mpc, one can calculate the SFCs of the 
mass field with Eq.(7). In these cases, Eq.(7) also insures that the variance of the density 
fluctuations 5 is less than 1, and therefore, the corresponding evolution of the mass field is 
still in the linear or quasi-linear regime. 



- 7- 



3. Lya forest estimator for rms bulk velocity 

3.1. Density-velocity relation in DWT representation 

By definition, a 1-D radial density perturbation S{r) is given by the projection of the 
3-D density perturbation 5{x,y,r) into 1-D, i.e. 

S{r) ^ I I d{x,y,r)dxdy = \ _ [ [ p{x,y,r)dxdy - p , (10) 

-t^x-t^y J La: JLy J-'xJ-'yP JLy 

where x Ly is the area of the coverage of samples on the celestial sphere. Both and 
Ly are larger than the scale considered in the radial (or redshift) direction. Since velocity 
field v(x) is statistically homogeneous, we have (l/LxLy) Jj^ {dvx{x,y,r)/dx)dxdy — 
0-1 Ly) /lJ^^(^2, y, r)-Vx{xi, y, r)\/ L^dy ~ 0. Similarly, {l/L^Ly) J^^ J^^{dvy/dy)dydx ~ 0. 
Thus, subjecting Eq.(l) to the operation (l/LxLy) ...dxdy, we have a 1-D equation 

where v{r) is a 1-D velocity defined by 

^^^^^JTlTl I '"r{x,y,r)dxdy (12) 
In redshift space, Eq.(ll) is 

where (1 -|- 2;) comes from the expansion of the universe, v is the radial component of the 

velocity field, and D = d.rjdz. For an LCDM model with + A = 1, we have D = 
c[HoE{z)]-\ H = HoE{z), E{z) = [Q(l + zf + A]^/^ (Peebles 1993), and / = f{z) « 
[0(1 + zfE{z)-^]^-^ (Lahav et al. 1991). 

The operator d/dz is almost diagonal in wavelet representation (Farge et al. 1996), and 
therefore, subjecting Eq.(13) to a DWT with scahng function (j)ji{z), we have 



1 , 2^-1 
1 + z 



HfD 



E4' / M^)^Mz)dz (14) 
i'=o 



where ej; is the SFC of the velocity field, given by ej; = J v{z)(j)ji{z)dz. The SFC ej; is 
proportional to the bulk velocity Vj on the spatial range given by the window function 4>ji, 
i.e. 

/ v{z)(f)ji{z)dz _ [¥ 
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Since J (j)ji(z)-^(f)ji'(z)dz — where is the so called connection coefficient of the 

basic scahng function (e.g. Restrepo & Leaf 1995), Eq.(14) can be rewritten as 



■' l'=0 



(16) 



This is the density- velocity linear Eq.(l) in the DWT representation. It is convenient for 



estimating e^i, from e^, and vice-versa. For D4 wavelet f^^i^ is non-zero only for | / — /'| < 2. 

For a Gaussian velocity field, we have (ejj/eji//) = if 1' 7^ 1" (Pando, Feng & Fang 2001), 
and therefore, Eq.(16) yields 



m\2\ 



(1 + ^)^ 



2^-1 



E 



=0 



(17) 



Moreover, for a uniform and isotropic random field, ((e^) ) and {{eji) ) are independent of 
position index Eq.(17) gives 

, -I 2 



m\2\ 



2^(1 + ^) 
LHfD 



(18) 



where W = Eu-r|<2(^^?i)'- 



3.2. Results of rms bulk velocity 

The ensemble of N Lya forests can be seen as a samphng of the sky with pencil beams. 
Because the operation (l/La-Ly) ...dxdy actually is an average over the celestial sphere, 

it can be statistically approximated by an ensemble average over the sky sampling. That is, 
the 1-D density perturbation 5{r) can statistically be estimated by 

1 1 ^ 

Thus, using Eqs.(7), (8), (9) and (12), we have 

f 1 ^ f 1 

1=1 

Substituting Eqs.(20) and (15) into Eq.(18), we finally have the estimator for the rms bulk 
velocity on scale j, crj'' = {vj)^^'^, of the underlying mass field of Lya forests 
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where R — cL/2'^ is the pixel size in km s~^, and the physical length scale correspond- 
ing to j is cL[2^ HoE{z)]~^ . The superscript Vr is to remind that it is a projection along 
the r-direction. The rms bulk velocity on a given scale is then determined from the flux 
fluctuations SF on the same scale. 

Using the estimator Eq.(21), we calculate the rms bulk velocities for the three redshift 
regions z ± 0.25 with z= 2.25, 2.75, and 3.25. The results are presented in Table 1. The 
scales in Table 1 have the same meaning as those in Fig. 1. In this calculation a is taken 
to be 1.6, and R = 293 km s^^. The factor A{z) is calculated from the mean transmission 
{F{z)) in each redshift range via A{z) = — lii{F{z)). 

It should be pointed out that the average over N pencil beams in eq.(19) will cause a 
Poisson error, which then contributes to the variance of S{r). This problem is the same as 
the Poisson correction in calculating the power spectrum of galaxies (Peebles 1980, Fang & 
Feng 2000). In Table 1 this correction has already been applied. 

The errors in Table 1 are calculated from different realizations of phase-randomized 
spectra which preserves the Gaussianity on large scales (Jamkhcdkar, Zhan, & Fang 2000; 
Zhan, Jamkhedkar, & Fang 2001). It is essentially similar to the bootstrap error estimation 
(Mo, Jing, & Borner 1992). Besides this error, there are also systematic errors caused by 
the uncertainties of the factor aA{z) in Eq.(7). For a given redshift range, the uncertainty of 
the mean transmission A{z) is less than 10%, and the uncertainty of a is about (1.9-1.5)/1.6 
= 25%. Therefore, the total systematic error would be about 30%. 

Continuum-fitting may cause error if the continua Fc of the 60 QSOs fluctuate signif- 
icantly on the same scale considered, say 92 h~^ Mpc. However, the PDFs of ej^ on large 
scales are perfectly Gaussian. Therefore, the SFG ej^ may be significantly contaminated by 
the continuum fluctuations only if the PDF of the continuum fluctuations is also Gaussian. 
This seems quite unlikely, as the emission from each QSO is highly non-thermal. 

All the 1-D rms bulk velocities listed in Table 1 are consistent with the assumption 
in §2, i.e. the peculiar velocity Vr is much smaller than the spatial size considered, i.e. 
Hor > 2300 km s~^. Moreover, the error from the resolution R of the samples is also small. 

Table 1: 1-D rms bulk velocities aj'~ (km s~^) 



Scales (h-^Mpc) 



92 46 23 



z = 2.25 
z = 2.75 
z = 3.25 



130±50 190±55 230±50 
62±40 86±30 94±20 
110±45 170±50 190±40 
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as R is also much less than the spatial size considered. In other words, the position and 
density uncertainties caused by Vr in Eq.(2) and R are negligible in comparison with the 
scales listed in Table 1. We should point out that the result of Table 1 is not sensitive to the 
choice of wavelets. Our analysis relies only on the orthogonality, completeness, and locality 
of the scaling function. Therefore, all wavelets with compactly supported basis will produce 
similar results. 

As expected. Table 1 shows that decreases with scales from 23 to 92 h"^Mpc, and it 
also decreases with redshift from 2.25 to 3.25. However, there is a dip of the bulk velocities 
on all scales at around z = 2.75. This should be verified with other sets of Lya QSO 
transmission fluxes. 



4. Discussions and conclusions 

To properly compare the estimated crj'' with other observed or theoretical results, we 
should emphasize that crj'' and Vj are defined by the DWT decomposition of velocity field 
(Eq. [15]), not by the velocities of particles or galaxies. It has been point out recently that 
the velocity defined by the particle-counting or galaxy-counting method may have different 
statistical properties from that defined by velocity field decomposition (Yang et al. 2001.) 
Therefore, we should compare our results with those that are also based on field decompo- 
sition if available. 

Prom Eqs.(12) and (15), we have 

_ II^ I "^rjx, y, z)(l)ji{z)dxdydz 
L^Ly j (j)ji{z)dz 

In terms of the DWT analysis, the average (l/Lx) Jj^ ...dx can be replaced by scahng func- 
tion projection window J ...(f)j^^i{x)dx/ J (f)j^^i{x)dx, where corresponds to scale L^- Thus, 
Eq.(22) can be rewritten as 

^ /^r-(^, y, z)(j)j^,i^{x)(pj^^i^{y)<f)ji{z)dxdydz ^ 

! 4>jaA^)K,iMU^)d^dydz ^ ^ 

Since the scales and Ly are much larger than the radial scale L/2-' considered, we have 
3x-,3y ^ 3- Therefore, actually is the 1-D rms bulk velocity of the 3-D mode {jx,jy, j) with 
jx-ijy ^ J, i-e. crj'' = {Vj^jy^y^'^ = '^j^,jy,j- For a homogeneous and isotropic random mass 
field with a relatively fiat power spectrum around 100-150 Mpc, we have a'"'', , ~ -s/Scr''' . . 
if jx,jy ^ j-i or equivalently, and Ly are much larger than the physical scales of j (Yang 
et al 2001). 



- 11 - 



Since the results listed in Table 1 are derived for high redshifts, there is no directly 
measured result available for comparison. The bulk velocity in a sphere of radius 50 h~^Mpc 
around the Local Group is found to be 370 ±110 km s~^ (Dekel et al 1998). The scale of 
the sphere is comparable to the linear scale 92 h~^Mpc in Tabic 1. If this result can be 
considered as the bulk velocity of a box 100^ h~'^Mpc^, the corresponding 1-D bulk velocity 
Vj on scale 100 h'^Mpc (with L^, Ly > lOOh^^Mpc) is about 370/^3/^3 = 123 km s^^ (the 
first y/S is due to the transform from 3-D to 1-D of the velocity, and the second y/S is due to 
the transform from mode {j,j,j) to (0, 0, j)). Considering hnear redshift evolution, this 1-D 
bulk velocity yields 99, 92 and 88 km s~^ at redshift z =2.25, 2.75, and 3.25, respectively. 
These results are consistent with that shown in Table 1. Nevertheless we should remember 
that the bulk velocity 370 km s~^ is from one realization, not rms value, and it may have 
the uncertainties of galaxy-counting as well as galaxy bias. 

The 1-D rms bulk velocity in the DWT modes (j,j,j) at redshift z — has been 
calculated from N body sinmlations. The results on 92 h~^Mpc are 60±10 km s~^, 70±20 
km s-\ and 80±20 km for the SCDM, LCDM and rCDM models respectively (Yang et 
al 2001). All models give lower rms bulk velocity than that shown in Table 1. One cannot, 
however, conclude that the N body simulated result is inconsistent with that shown in Table 
1. This is because the rms bulk velocity is sensitive to the perturbation on large scales, while 
the simulation box in Yang et al (2001) is 256^ h~^Mpc^. It has been pointed out (Tormen 
& Bertschinger 1996) that the hnear 3-D rms bulk velocity of a cube of side 100 Mpc is well 
over 500 km s~^ for an SCDM model with as — 1. Correspondingly, for an LCDM model 
with (78 = 0.7 the 1-D rms bulk velocities of DWT mode (L^;, Ly, 92 h^^Mpc), i.e. Co^oj' 
110 km s~^, 105 km s~^, and 100 km s"-*^ at z =2.25, 2.75, and 3.25 respectively. Therefore, 
the results shown in Table 1 basically are consistent with the linear evolution theory within 
la confidence level. 

We should also emphasize that the average over a QSO ensemble, i.e. average over 
the celestial sphere, is important. The errors in Table 1 also comes from the fact that the 
number of Lya forests in a given redshift bin is still small (large Poisson correction caused 
by eq. [19]). Nevertheless, one can already conclude that a set of Lya forests can provide 
a valuable estimate of the rms bulk velocity of the underlying mass field on scales from a 
few tens to < 100 h~^Mpc. With a large set of Lya forests, one may map the bulk velocity 
on large scales in a wide range of redshift, and then, set a robust constraint on the power 
spectrum on large scales. 

The authors would like to thank David Burstein for helpful suggestions. 
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Fig. 1. — Probability distribution functions of ej^s of the data set consisting of 
spectrum segments within the redshift range from 2.5 to 3.0. An LCDM {Q, — 
A = 0.7) universe is assumed to calculate the comoving length scales, ^ks is 
probability of the PDF being Gaussian. 
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